Imaging Pauli repulsion in scanning tunneling microscopy 
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A scanning tunneling microscope (STM) has been equipped with a nanoscale force sensor and 
signal transducer composed of a single D2 molecule that is confined in the STM junction. The 
uncalibrated sensor is used to obtain ultra-high geometric image resolution of a complex organic 
molecule adsorbed on a noble metal surface. By means of conductance-distance spectroscopy and 
corresponding density functional calculations the mechanism of the sensor /transducer is identified. 
It probes the short-range Pauli repulsion and converts this signal into variations of the junction 
conductance. 

PACS numbers: 68.37.Ef, 68.37.Ps, 67.63.Cd 



Since its invention the scanning tunneling microscope 
(STM) has become an important tool of nanoscience, be- 
cause it routinely provides Angstrom-scale image resolu- 
tion on various sample surfaces [1-6]. However, STM 
suffers from a serious drawback — the inability to resolve 
complex chemical structure. This disadvantage arises be- 
cause the STM probes the local density of states (LDOS) 
in the vicinity of the Fermi level, while details of the 
chemical structure are primarily encoded in lower-lying 
orbitals. Better access to chemical structure is therefore 
provided by mapping the total electron density (TED). 
Indeed, it has been shown recently that non-contact 
atomic force microscopy is able to resolve the inner struc- 
ture of a complex organic molecule, by imaging short- 
range repulsive interactions that depend on the TED 
[7]. Even earlier, however, it has been demonstrated 
that STM acquires very similar force imaging capabilities 
when operated in the so-called scanning tunneling hydro- 
gen microscopy (STHM) mode [8] — STHM images indeed 
closely resemble the chemical structure formulae of the 
investigated compounds (see ref. [8] and Fig. la). In this 
letter we present an analysis of the STHM junction by 
means of di/cH^z) -spectroscopy and density functional 
theory that allows us to explain its imaging mechanism 
for the first time. 

The experiments were performed on PTCDA/Au(lll) 
with a Createc low-temperature STM operated at 5- 
10 K in ultra-high vacuum (UHV). The electrochemi- 
cally etched W tips and the Au(lll) surface have been 
prepared using Ar + sputtering and temperature anneal- 
ing in UHV. The STM tips were additionally prepared 
by indentation into the clean gold surface. PTCDA 
molecules were deposited from a quartz crucible mounted 
in a home-built Knudsen cell. Deposition of H2 or D2 was 
performed according to the recipe described in ref. [8]. 
Since both H2 and D2 yield similar results, we restrict 
the discussion in this letter to D2. 

We start by summarizing the features of STHM that 
have been reported before [8]. The best STHM resolu- 




FIG. 1: (a) STM (top) and STHM image (bottom) of PTCDA 
(3,4,9, 10-perylenetetracarboxylic-dianhydride) on Au(lll): 
1.3x0.7 nm 2 , constant height, D 2 , V = 316 mV (STM) 
and V — —5 mV (STHM). The chemical structure formula 
of PTCDA is shown for comparison, (b) dl/dV spectrum 
measured in the center of PTCDA/Au(lll), recorded with 
lock- in detection (10 mV modulation, frequency 2.3 kHz). 
Go = 2jL_ = (12.9 /cQ) _1 is the quantum of conductance. 



tion (as shown in Fig. la) is achieved in constant-height 
mode with rather low tunneling bias \V\ < 10 mV. The 
appearance of STHM imaging after D2 dosing coincides 
with the emergence of non-linear differential conductance 
spectra G(V) = dI/dV(V) close to zero bias. In par- 
ticular, inelastic features appear (in Fig. lb at ±VJ ne i) 
[8-11]. In some cases, as in Fig. lb, a pronounced zero 
bias anomaly has also been observed [8]. Imaging in the 
presence of D2 can be switched reversibly between the 
STM and STHM modes by changing the applied bias V 
[8]: If \eV\ exceeds |eVi ne i|, the junction images in the 
conventional STM (or LDOS) mode. If, however, \eV\ 
is sufficiently smaller than |eVJ ne i|, the junction operates 
in the STHM mode, yielding images with ultra-high ge- 
ometric resolution. The nature of the inelastic features 
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has been investigated in detail in nano junctions contain- 
ing H2 [8-12]. They are assigned to transitions between 
two structurally different states of the junction. However, 
the precise nature of the associated two-level system is 
still debated [9-11]. 

Understanding the structure of the STHM junction 
is vital for identifying the imaging mechanism. Based 
on the conductance values in our experiments (> 5 x 
10 _5 Go) we conclude that the tip-surface separation 
must be smaller than 1 nm [13]. At such distances 
the junction can accommodate only a single monolayer 
of D2 [14, 15]. Clearly, the deuterium-induced imaging 
mode requires the presence of D2 just below the tunnel- 
ing tip apex. Given the high resolution that is achieved 
in STHM, which in Fig. la is of the order 50 pm, at most 
one D2 molecule can be in the active area of the junction. 
Therefore we can model the junction in the STHM mode, 
i.e. at \V\ < \Vine\\, by a single D2 molecule physisorbed 
[17] between the tip apex and the sample surface. For 
IV'I > I Vinei I ? the D2 molecule is displaced away from 
the tip apex, as proven by the recurrence of conventional 
STM imaging in conjunction with the occurrence of a 
structural transition at |Vi ne i|. 

Switching the bias between \V\ < \V[ ne i\ and \V\ > 
\V[ ne \\ we can study the properties of a given junction, 
and in particular the differential conductance G as a func- 
tion of tip-surface distance z, in presence and absence 
of the confined D2 molecule, without any other struc- 
tural changes of the junction. Fig. 2a displays the re- 
sults of G(z)-spectroscopy for the empty junction (curve 
measured with \V\ > |Mnei| and labelled G vac ) and for 
the junction with D2 (curve measured with \V\ < |Mnel| 
and labelled Gd 2 )- G vac increases exponentially with 
decreasing z, as expected for vacuum tunneling. In con- 
trast, Gd 2 behaves non-exponentially. Two opposing 
tendencies are observed: For intermediate z, Gjj 2 ex_ 
ceeds G vaci while for smaller tip-surface distances the 
situation reverses. Based on the data in Fig. 2a, we can 
define three characteristic regimes. 

At the shortest distances recorded in Fig. 2a (regime 
3) Gd 2 exhibits increased noise and conductance jumps. 
Such conductance changes usually occur when the tun- 
neling junction undergoes structural modifications due to 
the tip contacting the sample surface. Notably, the dis- 
continuities in Gd 2 occur at tip-surface distances where 
G vac still behaves strictly exponentially, i.e. the empty 
junction is still out of contact. Hence, the contact 
in question must occur via the confined D2 molecule. 
Clearly, the associated structural changes will occur in 
the softest part of the junction, i.e. the D2 molecule, 
which eventually is squeezed out of the junction. Indeed, 
images measured at the onset of regime 3 (images 5 and 
6 in Fig. 2c) show a sudden loss of STHM resolution. 

Having identified regime 3 with the squeezing of the 
D 2 molecule out of the junction, we can associate the 
preceding regime 2 with the gradual compression of the 




FIG. 2: (a) ^7(2) spectra recorded at fixed bias, measured 
with D2 on PTCDA/Au(lll) (z-axis scale given by scale bar). 
Go — ^f- is the quantum of conductance. Gd 2 (magenta): 

V = -5 mV. G vac (black): V = -130 mV. All spectra 
recorded from the same stabilization point, (b) Rg(z) curves 
(cf. text) measured at positions r\ (blue) and r<z (red) marked 
in the lower inset. The ratio G vac (ri) /G vac (r2) is shown in 
green, (c) STHM images (1.3 x 1.3 nm 2 , constant height, 

V = -5 mV) of PTCDA/Au(lll) measured with D 2 at dif- 
ferent z indicated by black triangles in panel b. 



junction that eventually causes this squeeze-out. To 
quantify the effect of D2 on the junction conductance 
in regimes 1 and 2, we define the conductance ratio 
Rq(z) = Gd 2 (z) / Gv&c(z) • As can be seen in Fig. 2a, 
Rq(z) can be larger or smaller than 1, depending on 
the value of z. Fig. 2b (red, blue) shows that in regime 
2, where the best STHM resolution is recorded (images 
3 to 5 in Fig. 2c), Rg(z) curves measured at different 
lateral positions above the PTCDA molecule each dis- 
play a distinct slope. At the same time G vac curves do 
not vary appreciably from point to point above PTCDA 
(Fig. 2b, green), which is consistent with the blurred and 
featureless STM images that are recorded with the empty 
junction at \V\ > |Vmei| (Fig. la). The Rq(z) curves in 
Fig. 2b hence show that the STHM contrast in Fig. la 
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FIG. 3: DFT-LDA simulated n t ,D 2 (E F , z') /n t , vac (E F ) vs. 
D2-tip distance z for s- (black) and p-type (red) orbitals at 
the Au atom below the deuterium molecule. n t ,D 2 (^F, z f ) is 
the LDOS of the model tip at the Fermi level at given D2-tip 
distance z' and n t)V ac(^F) is LDOS at the Fermi level of the 
bare tip electrode. Equilibrium distance z' e(l is indicated by 
the vertical line. For details of the simulation, cf. ref. [18]. 



can be ascribed to the effect of D2 on the junction con- 
ductance which becomes pronounced in a narrow range 
of z -values in regime 2. We therefore need to study the 
nature of regime 2 in more detail. 

To investigate the influence of D2 on the tunnelling 
conductance of an STHM junction in regime 2, we have 
carried out a density functional theory (DFT) calcula- 
tion in which we systematically varied the distance z' 
between D2 and an Au(lll) surface (Fig. 3 inset) [18]. 
This surface was chosen to model the tip electrode in our 
STHM experiments. We calculate the model tip DOS 
at the Fermi level, ntp 2 (Ep, z'), as the LDOS at the 
metal atom located directly underneath the D2. Divid- 
ing nt,D 2 (^) z') by n t ,vac(^F), the LDOS of the bare 
tip without D 2 , we find that the tip DOS decays sub- 
stantially with decreasing z' (Fig. 3). The origin of this 
behavior is the Pauli exclusion principle: To minimize 
overlap between the closed shell of D2 and the metal elec- 
trons, both of their wave functions must rearrange locally, 
which depletes the metal's local DOS in the vicinity of 
the Fermi level [15, 19, 20], while the associated energy 
cost leads to a repulsive force between D2 and the metal 
(Pauli repulsion). 

In the limit of low tunneling bias TersofT-Hamann the- 
ory of STM predicts G oc n t ^ 2 (E Fl z / )n s (E Fl f t ) 1 where 
n s is the sample LDOS at tip position f t [21]. Ac- 
cordingly, G and Rq must decrease proportionally to 
nt^ 2 {EF-,z') as D2 approaches the tip. DFT results 
shown in Fig. 3 suggest that the rate of tip DOS decrease 
should be in the range from 0.2 to 1 A .At the same 
time experimental data from Fig. 2b show rates between 
1.1 and 1.5 A -1 . To be able to compare both results 
we additionally have to divide experimental values by 
the factor ^- accounting for the differences between the 



scales z and z' . must be in the range between ~ 
(soft sample) to 1 (hard wall sample). Given the effects of 
the unknown tip shape and the neglect of modifications 
of the sample LDOS by Pauli repulsion [22], the agree- 
ment between experiment and simulation is remarkable 
and we can thus conclude that the z-variation of Rq{z) 
in regime 2 can be explained as an effect of Pauli repul- 
sion in the STHM junction. 

This result holds the key for understanding the con- 
trast formation in STHM. To demonstrate this, we first 
discuss the STHM contrast above an inherently simple 
object, namely Au adatoms on Au(lll), before turning 
to the more complex PTCDA molecule with its inter- 
nal structure. Fig. 4 (bottom panel) shows the experi- 
mental image of an Au adatom dimer. The image has 
been recorded at constant height, nearly zero bias volt- 
age (2 mV), and with a junction containing D2, i.e. under 
STHM conditions. We observe two well-separated struc- 
tures, each of which corresponds to one of the adatoms. 
In comparison to the flat sample surface, the adatoms ap- 
pear bright, i.e. with a large G(r t ), because the sample 
LDOS n s (~ Ef, f t ) at tip positions close to the adatoms, 
e.g. at f t2 , is increased with respect to the one at r t i, due 
to a reduced effective tip-sample distance (cf. the top 
panel of Fig. 4). So far this is not different from conven- 
tional constant height imaging in STM. However, in the 
center of each of the adatoms (i.e. at r t ~ fts) a dark area 
is observed in Fig. 4. In contradiction to conventional 
STM [23], Fig. 4 clearly implies that G(r t3 ) < G(r t2 ). 
The reason for this deviation from the normal STM be- 
havior can be found in the presence of D2 in the junc- 
tion, and in particular in its trajectory, which is displayed 
schematically in Fig. 4 (top): As the tip moves from ft 2 
to ft3, the D2 molecule will have to move to a new vertical 
equilibrium position closer to the tip (smaller z') because 
at that position the stronger Pauli repulsion from the 
adatom is balanced by a stronger Pauli repulsion from 
the tip. In conjunction with Fig. 3, this must lead to a 
sharp reduction in nt(— E F ). According to Fig. 4 this 
reduction overcompensates the rise in n s (^ Ep,ft) [22] 
and leads to the dark areas in the centers of the adatoms. 
The analysis of the adatom image thus shows that the 
STHM contrast can be understood as an (x,y)-map of 
the short-range Pauli repulsion from the sample surface 
acting on the D2 molecule in the STHM junction, super- 
imposed over the conventional LDOS contrast. 

With this knowledge, we can finally analyze the STHM 
contrast generation above PTCDA. As in the case of the 
dimer, the D2 molecule follows the tip and probes lat- 
eral variations of the Pauli repulsion from the adsorbed 
PTCDA. For example, when the tip moves from a posi- 
tion above the center of a C6 ring to a position directly 
above a carbon atom, the D2 molecule in the junction 
will — similarly to the trajectory shown in Fig. 4 — move 
to a higher equilibrium position closer to the tip, be- 
cause of the increased TED above the carbon atom. As 
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FIG. 4: (a)(b) STHM image of a Au dimer on Au(lll) (con- 
stant height, V —2 mV, D2 ) (bottom panel) and schematic 
sketch of contrast generation. Cf. text for details. 



in the case of the adatoms discussed above, this leads 
to a reduced n t (~ Ep) and therefore lower conductance. 
In the STHM image the carbon atoms of PTCDA (and 
by a similar argument the a-bonds between the carbons) 
therefore appear darker than the "empty" spaces inside 
the C6 and C5O rings of the PTCDA backbone, just as 
observed in the image of Fig. la. 

In conclusion, we arrive at the following model of 
STHM imaging: A single D2 molecule is physisorbed in 
the STM junction, such that it is confined directly un- 
derneath the tip apex. This molecule is the crucial el- 
ement in the STHM imaging process, as it probes the 
short-range Pauli repulsion from the surface (sensor ac- 
tion) and transforms this force signal into variations of 
the tunneling conductance (transducer action) , the latter 
again via Pauli repulsion. Because of its nanoscale size, 
the sensor is insensitive to long-range forces. Clearly, the 
sensor /transducer modulates the tunneling conductance 
on top of the normal LDOS contrast. As long as the 
Pauli-induced conductance modulation is larger than the 
LDOS-induced change of the background conductance 
itself, the image will be dominated by the STHM con- 
trast. We note that the described functionality should 
also work with other closed-shell particles besides hy- 
drogen and deuterium [24]. Comparing STHM to con- 
ventional STM, direct tunneling between tip and sample 
surface still forms the basis of imaging in STHM. How- 
ever, in STHM a compliant element that is sensitive to a 
laterally varying sample property other than the LDOS 
is added to the tunneling junction, where it modulates 
the tunneling current that is used for imaging. 
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